Three-body systems with attractive 1/r potentials 
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We have used the hyperspherical adiabatic representation to describe the system of three identical 
bosons in an spin stretched state interacting by an attractive 1/r potential. A proposal has been 
made how such a system might be realized experimentally in cold trapped atoms using extremely 
off-resonant laser fields [Phys. Rev. Lett. 84, 5687 (2000)]. We have obtained effective potentials, 
channel functions, and nonadiabatic couplings for this gravity-like interaction, allowing us to calcu- 
late the ground state energy with accuracy that substantially improves upon previous results. We 
have similarly calculated the energies for the first four + excited states. These results show that 
the simple adiabatic hyperspherical approximation offers an accurate description for such a system. 

PACS numbers: 31.15.Ja,32.10.-f,34.20.Cf 



I. INTRODUCTION 

Recently, a scheme for inducing gravitation-like in- 
teratomic potentials has been proposed 0, opening up 
the possibility to create self-bound Bose-Einstein con- 
densates (BECs) H S S H H. In such a scheme, 
the gravitation-like interatomic potential can be achieved 
by irradiating the atoms with intense, extremely off- 
resonant electromagnetic fields. The usual strong 
anisotropy due to dipole-dipole interactions can, in fact, 
be averaged out Q by the proper combination of laser 
beams, leaving a —ujr potential, where u is the strength 
of the potential (the analogue of GMm, where G is the 
Newton's constant and M and m are the masses) and r 
is the interparticlc distance. The strength u of the po- 
tential can be adjusted by changing the laser intensity 

Q 

In ultracold atomic gases, two interesting regimes for 
self-bound BECs have been predicted, assuming that 
the short-range interatomic interactions can be indepen- 
dently tuned by, say, applying a magnetic field near a 
Feshbach resonance In one regime, the attractive 1/r 
interactions are balanced by the repulsive mean field in- 
teration assuming a positive two-body scattering length 
and negligible kinetic energy In the other regime, the 
balancing factor is kinetic energy, assuming negligible 
mean field interactions. In both regimes, the resulting 
BEC is self-bound. From a broader point of view, the 
induced gravitation-like interaction might make possi- 
ble experimental emulation of boson stars (a system of 
self-gravitating bosons) in the regime where the kinetic 
energy balances the — ujr potential 0, ^j. Moreover, 
purely attractive 1/r potentials constitute an interesting 
contrast to the attractive and repulsive Coulomb poten- 
tials atomic physics are used to. 

The existence of a lower bound for the ground state 
energy in many-body systems interacting through attrac- 
tive 1/r potentials is of fundamental importance in or- 
der to prove the existence of the thermodynamical limit 
and the stability of normal matter |ll| . It was shown 
in Ref. 0] that for a system of N identical, spinless (or 
spin stretched) bosons of mass m interacting gravitation- 



ally, the lower and upper bounds for the ground state 
energy are, respectively, — jqN 2 (N — 1) G 2 m 5 /h 2 and 
-0.05427V(7V- 1) 2 G 2 m 5 /h 2 , where the upper bound was 
obtained variationally. For small N, however, the dis- 
crepancy between the lower and upper bounds becomes 
large. For TV — 3, using a more refined trial function, 
they obtained -0.95492 G 2 m 5 /H 2 for the upper bound, 
representing a difference of about 15% between the up- 
per and lower (— | G 2 m 5 /h 2 ), bounds and a ground state 
energy equal to E = —1.067 G 2 m 5 /H 2 . 

In this paper, we have used the adiabatic hyperspheri- 
cal representation for this systm to obtain effective three- 
body potentials, the corresponding channel functions, 
and the nonadiabatic couplings. Using these, we calcu- 
late the ground state and low-lying + excited state en- 
ergies converged to seven digits. We have also used the 
adiabatic hyperspherical representation to obtain lower 
and upper bounds that differ by about 0.1%, indicating 
that a simple single-channel description offers a quite ac- 
curate description of such systems. 



II. THE ADIABATIC HYPERSPHERICAL 
REPRESENTATION 

We have solved the Schrodinger equation in hyper- 
spherical coordinates. After separation of the center-of- 
mass motion, the system is described by the hyperradius 
R which gives the overall size; three Euler angles a, (3 
and 7, specifying the orientation of the plane containing 
the three particles relative to the space-fixed frame; and 
other two hyperangles if and 9, describing the internal 
relative motion between the particles. We have defined 
if and 9 as a modification of Smith- Whitten coordinates 
0,0,0]. The key to the adiabatic hyperspherical rep- 
resentation is that the dynamics of the three-body system 
is reduced to collective motion under the influence of one- 
dimensional effective potentials in R, which is governed 
by a system of ordinary differential equations. 

The hyperspherical coordinates are introduced through 
the mass-scaled Jacobi coordinates p\ and pi (see Fig. ^) 
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FIG. 1: The mass-scaled Jacobi coordinates for systems with 
three particles. 



defined as 



Pi = (72 - n)/d, 

mi?x + m 2 r 2 



92 



'"3 



mi + m 2 



(1) 



In the above equations is the position of the particle 
i (of mass mi) relative to a space-fixed frame. For three 
identical particles of mass m, we define a three-body re- 
duced mass as /i = m/V3 which gives d = 2 1 / 2 /3 1 / 4 
[T3. ITHI|. It is important to note that the hyperradius, 



R 2 = P i 



Pi 



Re [0,00), 



(2) 



is an invariant quantity, i.e., it does not depend on the 
particular choice of the hyperangles or labels of the par- 
ticles. 

The Schrodinger equation can be more conveniently 
written in terms of the rescaled wave function ip — 
i? 5 / 2 <F, as 



h 2 d 2 

2^ OR 2 



Had(R, SI) 



i>(R,Sl) = E1>{R,fl), (3) 



where E is the total energy and H ad is the adiabatic 
Hamiltonian given by 



H ad (R,n) = 



2[iR 2 



a 2 (si) + ^ 



+ V(R,<p,9). (4) 



The adiabatic Hamiltonian H ad contains all hyper- 
angular dependence, represented collectively by SI = 
{(p, 9, a, f3, 7}, and includes the hyperangular kinetic en- 
ergy in the grand angular momentum operator A 2 as well 
as all interparticle interactions V. 

In the adiabatic hyperspherical representation, the to- 
tal wave function is expanded in terms of the channel 
functions <I>„(i?;fi), 



4> n (R,tt) =J2F nu {R)$„(R;n), 



(5) 



where F nu (R) are the hyperradial wavefunctions, n labels 
the different energy eigenstates for a given v, and v repre- 
sents all remaining quantum numbers necessary to spec- 
ify each channel. The channel functions <!>„(-/?; SI) form a 
complete set of orthonormal functions at each value of R 
and are eigenfunctions of the adiabatic Hamiltonian: 



H ad (R, n)$ u (R; Si) = U V {R)® V (R; SI). 



(6) 



The eigenvalues U V (R) help define effective three-body 
potentials for the hyperradial motion. 

Substituting Eq. (JjjJ into the Schrodinger equation 
(|3"|) and projecting out $„', we obtain the hyperradial 
Schrodinger equation 



h 2 d 2 
T^idR 2 



U V {R) 



2P UV .{R) 



dR 



Qw'{R) 



F V .{R)=EF„(R), 



(7) 



r 



which describes the motion of the three-body system 
under the influence of the effective potentials U V {R) — 
Qw{R)I^P- The nonadiabatic coupling terms P VV >{R) 
and Q VV ' (R) drive inelastic collisions three-body scatter- 
ing processes and are defined as 



P vv r{R) = (($ 



d 




dR 





and 



Q VV ,{R) = (($ 



dR 2 



(8) 



(9) 



where the double brackets denote integration over the 
angular coordinates SI only. As it stands, Eq. J7J is exact. 



In practice, of course, the sum over channels must be 
truncated. In fact, the accuracy of the solutions can be 
monitored with successively larger truncations since the 
bound state energies obtained at each stage are an upper 
bound by the variational principle. 

In this paper, we explore the solutions of the system of 
differential equations J7J for three particles with attrac- 
tive X/r interactions. We determine the effective poten- 
tials and couplings by solving Eq. |JBJ for the J* = + 
symmetry, where J is the total orbital angular momen- 
tum and 7T is the total parity. The low-lying bound state 
energies are then determined by solving Eq. J7J. 

In order to solve the adiabatic equation |JBJ, we have 
expanded the channel functions ^^(i?; SI) in terms of the 
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Wigner D functions 0,13113, 

$f l7r (R; Q)=J2 ^Ku(R; 0, tp)D J KM (a, 0, 7 ), (10) 

K 

where K and M are the projection of J onto the body- 
fixed and space- fixed z— axes, respectively. After pro- 
jecting out the D functions, the resulting coupled system 
of partial differential equations for <pKv{R] 9, tp) is solved 
(for each value of R) by expanding 4>Ku(R\ 9, tp) on a 
direct product of fifth order basis splines 18] in the hy- 
perangles 9 and tp 0,0]. For J 7r = + , of course, the 
sum involves only one term, requiring the solution of a 
single two-dimensional partial differential equation. 

The potential V in Eq. is given by a pairwise sum 
of attractive 1/r potentials, 

V{R,6,<P) = -—- — -—, (11) 
ri2 r 23 r 3 i 

where u is the gravitation-like coupling. The interparti- 
cle distances are given in terms of the hyperspherical 
coordinates by 

ri2 = 3" 1/4 i?[l + sin6»sin(^-7r/6)] 1/2 , 

r 23 = 3- 1/4 R[l + sm9sm(<p-5Tr/6)] 1/2 , 

r 3 i = 3" 1/4 i?[l + sin0sin(^ + 7r/2)] 1/2 . (12) 

Figure [3 shows the potential V(R,9,ip) as a function 
of 9 and tp at R = 100. The singular points at 9 = tt/2 
and tp = 7r/3, 7r and 57r/3 are the points where r 23 = 0, 
r 31 = and r 12 = 0, respectively. Notice that our choice 
for the hyperangles 9 and tp [14. llSl ] — like any of the so- 
caller democratic or Smith- Whitten coordinates — im- 
plies a periodicity of the potential in ip for identical par- 
ticles, which substantially simplifies the numerical solu- 
tion of Eq. © when symmetrizing the wave function. 
The two-dimensional equation for J 7 ' = + thus needs 
to be solved only from ip = to 27r/3, with the require- 
ment that the derivative of 4>(R; 9, tp) with respect to tp 
is zero at each boundary, followed by postsymmetriza- 
tion to extract the completely symmetric solutions. We 
have treated the cusp due the 1/r divergence as r — > 
by including more spline functions 0] at tp = ir/3. For 
other J v the solutions can be determined in a similar way 
00 

In the atomic case, suggested in Ref. [l|, u is given by 
(in S.I. units) 



v(R,e,<?) 



11 Iq 2 a 2 p 
4tt ce 



(13) 



where / is the laser intensity, q the photon wave number, 
and a p the atomic dynamic polarizability. So, under the 
conditions discussed in Ref. pj, the strength of the in- 
teraction is controllable via the laser's parameters. For 
the present calculations, however, we base our units on u, 
thus producing a unitless equation. Our length units are 
2S 2 /mu, and our energy units are mu 2 /2h 2 . These yield, 
in analogy to atomic units, a two-body energy spectrum 
En = -l/2n 2 . 
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FIG. 2: The potential V(R, 0, <p) at R = 100. Due to the 
symmetry properties of the three identical particles, the adi- 
abatic equation [Eq. @] is solved only for 6 = to 7r/2 and 
tp = to 2tt/3 



III. RESULTS AND DISCUSSION 

Figure 01 shows the effective three-body potentials in 
the form [— 2Z7„(i?)] -1 / 2 such that for large values of R 
they converge to the principal quantum number n 2 b asso- 
ciated with the two-body hydrogen-like subsystems (for 
R — > oo one particle is far from the others). The low- 
est potential in Fig. El converging to n 2 ^ = 1, supports 
the + bound states (or S e , in analogy with atomic spec- 
troscopic notation). In particular, it contains both the 
ground state and the S e series of singly excited states, 
using the language of atomic structure. The higher po- 
tentials (converging to n 2 \> > 1) support series of dou- 
bly excited states that are coupled to and can decay 
to the continuum of the lowest potential and are thus 
metastable. These comments draw on the similarity of 
these potentials to those for three-body systems like He 
and H _ [19J. There are minor differences, of course, due 
to the different permutational symmetries of the two sys- 
tems and to the absence of the Coulomb repulsion for the 
present case. 

Since the J 77 = + channel functions depend only on 
the hyperangles 9 and tp they can be plotted in their 
entirety for each value of R. Figures 0] and [5] show the 
channel functions for the two lowest channels at R = 0.69 
and 100 for v = 1, and R = 5.75 and 100 for v = 2. The 
first R value lies near the respective potential minima; 
and the second, in the asymptotic region. For small R, 
we plot <!>„ in the range < tp < 2ir/3 from which the 
function in the whole range < tp < 2tt can be obtained 
by symmetry (translation of the plotted portion by 27r/3 
and 47r/3). For large R, however, we plotted the solution 
only in the range 7r/6 < tp < tt/2 to emphasize the two- 
body character of the solution. 

The hyperangular distributions are useful because they 
reveal the geometry of the system. At 9 = 0, for instance, 
the atoms form an equilateral triangle, while for 9 = n/2 
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FIG. 3: The hyperspherical potentials for three identical 
bosons with attractive 1/r interactions. For large values of 
R, [— 2U V (R)]~ 1 ^ 2 converges to the principal quantum num- 
ber n2b for the hydrogen-like subsystems. 



they lie along a line. Figure EI a) shows that near the 
v = 1 potential minimum (see Fig. |2J) the lowest channel 
function is spread out over the entire hyperangular plane 
with an increased amplitude near the two-body coales- 
cence point (r 3 i = 0) at = n/2 and ip = n/3. This 
point corresponds to a linear configuration with two of 
the particles closer to each other than to the third par- 
ticle, representing strong two-body correlations. Since 
this R is the equilibrium distance, the three particles in 
the ground state thus assume all triangular shapes be- 
tween linear and equilateral with a preference for linear 
shapes having two particles close to each other (taking 
into account the sin29 volume element). As R increases 
[Fig.Qfb)] the channel function "collapses" to the region 
around the coalescence points, displaying the two-body 
character it must have for R — > oo — in this case the 
Is state. Figures Efa) andOJb) show the v = 2 channel 
function. They show much the same behavior as v = 1 
except with the necessary addition of a node. Figure 
Etb) shows that v — 2 converges to the 2s state of the 
two-body subsystem. Note that for identical spinless (or 
spin-stretched) bosons, 2p two-body states are not al- 
lowed by symmetry, so there is only a single potential 
correlating to ri2b = 2 at R — > oo. 

In Figure we show the nonadiabatic couplings 
Pi/v'(R) between the lowest channel and the next three 
channels [P vv = and P vv i = —P v i u ) calculated from 
Eq. ©. The fact that the coupling Px% in Fig. El is sub- 
stantially larger than the couplings with higher channels 
implies rapid convergence for the bound state energies as 
a function of the number of channels. Further, P\i peaks 
around R — 6 which correlates roughly with the loca- 
tion of an avoided crossing between the corresponding 
potential curves as expected (see FigEJl - 
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FIG. 4: The lowest J n = + channel function (y = 1) as 
function of 9 and ip at (a) R = 0.69 and (b) R = 100. 



We have solved Eq. Q and determined the bound state 
energies and hyperradial wavefunctions F nv (R) includ- 
ing up to 15 channels. For example, Figure [7] shows the 
ground state wave function for a three channel calcula- 
tion. The v = 1 component is associated with the lowest 
adiabatic channel, while the v = 2 and v = 3 are re- 
lated to the next two adiabatic channels and are present 
due to the coupling between the channels. Note that the 
probability given by V w = J °° \F nv {R)\ 2 dR due to the 
first term dominates the other contributions. In Tabled 
we show V nv= \ for the ground state (n — 1) and the 
four lowest excited states for a 15 channel calculation. 
The v = 1 adiabatic channel represents roughly 99% or 
more of the probability for each state, showing that the 
adiabatic expansion is, in fact, quite good. 

In Table UTI we show the ground state energy as a func- 
tion of the number of channels, further demonstrating the 
expected rapid convergence of the adiabatic expansion. 
The ground state energy for this system has certanly been 
calculated before. In Ref. 12], for instance, a ground 
state energy of E = -1.067 G 2 m 5 /h 2 = -2.134 mu 2 /h 2 
(in our units) was obtained. Comparison with Table [n] 
shows that since both calculations are variational, our 
single channel calculation already gives a more precise 
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FIG. 5: The first excited J w = 0+ channel function (y = 2) 
as a function of 9 and <p at (a) R = 5.75 and (b) R — 100. 




FIG. 6: Nonadiabatic coupling P vu i between the lowest chan- 
nel and the next three channels. 



result. We speculate that the large differences in the 
potential energies shown in Fig. [3] are the main reason 
that a single channel already gives such a good result. 
In fact, this channel separation is closely related to the 
small magnitude of the coupling terms shown in Fig. 
Table [H] also shows that our six channel approximation 
for the ground state gives a result converged to seven 




FIG. 7: First three components of the ground state hyperra- 
dial wavefunction F w {n — 1) associated with the first three 
channels v = 1, 2, and 3. 



TABLE I: The probability Vnv associated with the v = 1 
term of the expansion JSJ for the ground state (n — 1) and 
the next four excited states. 



n 


Pnl=/ n " , |Ji , nl| 2 dB 


1 


99.98% 


2 


99.80% 


3 


99.15% 


4 


98.90% 


5 


99.21% 


digits. 




TABLE II: Convergence of the ground state energy as a func- 


tion of the number of channels included in Eq. J7J . 


Number of channels 


Ground state energy ( mu 2 /h 2 ) 


I 


-2.136 033 


2 


-2.136 481 


3 


-2.136 523 


4 


-2.136 525 


5 


-2.136 526 


6 


-2.136 527 


15 


-2.136 527 



In Tabic II I II we give our converged results for the 
ground and first four excited states (calculated with 15 
channels). Is well known [2(j that the hyperspherical 
energy obtained disregarding all couplings and the en- 
ergy obtained considering only the diagonal coupling 
in Eq. (JJJ) are, in fact, lower and upper bounds for 
the ground state energy, respectively. In these ap- 
proximations, we obtained a lower bound correspond- 
ing to —2.138 650 mu 2 /h 2 and an upper bound of 
— 2.136033mu 2 /fi. 2 . The difference between them is about 
0.1% while the results obtained in Ref. 01 gi ye a differ- 



6 



ence of about 10%. By comparison, for a system like 
the He atom , where the electronic repulsion plays an 
important rule, the relative difference between the lower 
and upper bounds estimated from hyperspherical poten- 
tial curves is about 1%. 

TABLE III: Ground state and excited states energies E n ,v 
[y = 1) calculated using 15 coupled channels. 



n E niV ( mu 2 /h 2 ) 

1 -2.136 527 

2 -1.145 881 

3 -0.786 454 

4 -0.661 162 

5 -0.603 740 



IV. SUMMARY 

We have used the adiabatic hyperspherical represen- 
tation to describe system of three identical bosons with 
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